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Abstract. This article examines dynamical systems on a class of K3 surfaces in X with an infinite 
automorphism group. In particular, this article develops an algorithm to find Q-rational periodic points 
using information modulo p for various primes p. The algorithm is applied to exhibit K3 surfaces with 
Q-rational periodic points of primitive period 1, . . . , 16. A portion of the algorithm is then used to 
determine the Riemann zeta function modulo 3 of a particular K3 surface and find a family of K3 
surfaces with Picard number two. 



1. Introduction 

Let : F C -^V he a. morphism on a variety V . We denote the n*^ iterate of <j) as 

We say that P is preperiodic if {0'(P) | i £ N} is finite. In particular, we say that P is a periodic 
point of period n for (j) if 0"(P) — P and of primitive period n if, in addition, for all m < n, (j)™{P) ^ P. 

By work of Northcott there are only finitely many rational periodic points for a morphism of 
projective space. In particular, this implies that the possible primitive period of a rational periodic point 
is bounded. Morton and Silverman [10] conjecture for morphisms of projective space that this bound 
should depend only the degree of the number field, the dimension of the space, and the degree of the 
morphism. Fakhruddin [7] has shown that this conjecture implies a similar statement about morphisms 
of certain algebraic varieties. This article examines dynamical systems on a class of surfaces in x 
with an infinite automorphism group studied by Wehler [T5|. Silverman [131 Theorem 1.1] shows that 
the set of rational (pre)periodic points for the composition of the two non-commuting involutions on 
these K3 surfaces is a set of bounded height and hence finite. Call and Silverman [51 Appendix] give 
explicit formulas for computing the two non-commuting involutions on these K3 surfaces. This article 
addresses the following topics: 

(1) find Q-rational periodic points on Wehler K3 surfaces using information modulo p for various 
primes p G Q. 

(2) determine the number of rational points on a Wehler K3 surface in the field with elements, 
denoted ¥pm . 

In the course of addressing these topics, we produce an explicit family of K3 surfaces with Picard number 
two. Note that van Luijk [14] constructs explicit K3 surfaces with Picard number one and infinitely many 
rational points using more general methods than discussed in this article. However, since his examples 
have an elliptic fibration his methods are not wholly applicable to Wehler K3 surfaces. Additionally, the 
orbits of curves and points have been studied on these K3 surfaces by Baragar in [1] [3] . 

This article proceeds as follows. In Section [21 we define and recall some basic information about 
Wehler's class of K3 surfaces. In Section [31 we address finding Q-rational periodic points from data 
modulo primes. The key fact used to find rational periodic points is that if there is a rational periodic 
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point, then the primitive period of its reduction modulo a prime of good reduction must divide the 
primitive period of the rational point. While this method may return false positives due to checking 
only a finite number of primes, it provides an effective way to search for rational periodic points. This 
algorithm is used successfully in Pari/gp to find surfaces with Q-rational points of primitive period 
1, . . . , 16. In Section H we use the number of Fpm rational points with m = 1, . . . , 11 for a particular 
Wehler K3 surface to determine the Riemann zeta function of the surface using the Weil conjectures. 
An example is computed for p — S and shown to have Picard number two which produces an infinite 
family of Wehler K3 surfaces with Picard number two. Much of this work is from the author's doctoral 
thesis [51 Chapters 6 and 7]. 

There are also several similarly constructed classes of surfaces in Baragar [2] , Billard [4] , and McMullen 
[9] that may allow a similar type of explicit investigation. 

A PARI/gp script to search for a Q-rational primitive periodic point on a Wehler K3 surface and to 
check for the possible primitive periods of a periodic point on a Wehler K3 surface is available on the 
author's homepage. 



2. Wehler's K3 surfaces 

A Wehler K3 surface 5 C x is a smooth surface given by the intersection of an effective divisor of 
degree (1,1) and an effective divisor of degree (2,2). In other words, let {[xo,xi,X2], [yo,yi,y2]) = (x, y) 
be the coordinates for x P^, then S is the locus described by L(x, y) — Q{x,y) — for 

L{^,y) = X! ^v^iVj Q{^,y)^ X! bijkiXtXjykyi. 

0<i.j<2 0<i,jMd<2 

Theorem 1. |15[ Theorem 2.9] A general K3 surface formed as the vanishing locus of a degree (1, 1) 
and a degree (2, 2) effective divisor has Picard number two and an infinite automorphism group. 

When working with a particular surface, we will specify the vectors 

L = [aij] = [aoo,aoi,ao2,aio, . . . ,022] and 

Q ~ [bijkl] ~ [&OOOO1 ^OOOli ^0002, &0010> ■ ■ ■ 7 &2222]- 

The natural projections 

:P2xp2^p2, py:FlxPl^Fl 

induce two projection maps 

p,:S^Fl, py-.S^Fl. 

The projections p^ and Py are in general double covers, allowing us to define two involutions of 5, ax 
and ay, respectively. To define a.^, we consider p^^(x) = {(x, y), (x,y')} and define tTa;((x,y)) = (x, y'). 
To define ay, we consider Py^{y) = {(x, y), (x', y)} and define cry((x, y)) = (x', y). The maps a.^ and ay 
are in general just rational maps. However, if S* = V{L, Q) is smooth. Call and Silverman [21 Proposition 
1.2] show that ax and ay are morphisms of S if and only if S has no degenerate fibers, fibers of positive 
dimension. We call a surface with no degenerate fibers a non- degenerate surface. 

Remark. In the case where S is smooth and degenerate, Baragar [3J Section 2] shows how to extend ax 
and ay to morphisms defined also on the degenerate fiber. These extended morphisms are not treated 
here, but offer an interesting area of future study. 
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We adopt the notation from |5J and define 

L'j = the coefRcient of yj in L(x, y), 

= the coefficient of xj in ij(x, y), 

Qli = the coefficient of y^yi in (3(x, y), 

Q^j = the coefficient of XiXj in (5(x, y), 

-"ij — ^^i^j^kk ^ ^i^k^jk ^ ^j^kWik + l-'^fcJ Vij- 

for (i, J, fc) some permutation of the indices {0, f , 2} and * replaced by either x or y. 

3. Periodic Points 

We use good reduction information to develop an algorithm to find Wehler K3 surfaces with Q-rational 
periodic points of specified primitive period. 

Theorem 2. Let K he a number field. Assume that 5 C x is a Wehler K3 surface defined over 
K . Let (j) ~ ax ° <yy G Aut(S'), let P G S{K) be a periodic point for (j), and let (Z K be the set of 
primes where (p has good reduction. Define 

n — primitive period of P for (j), and 

mp — primitive period of P for (f> for p G *p. 

Then rup \ n for all primes p G *P and n = mp for all but finitely many primes p G 
Proof. To show divisibility we know for primes of good reduction 

0"(P) ^TiP) =P 

so we have rup < n. Since the smallest period is the greatest common divisor of all of the periods, we 
must have mp \ n. 

For the second statement, let Xi{P) and yi{P) for < i < 2 be the coordinates of P. The assumption 
that (P) = P is equivalent to the system of equations 

CiXo(<^"(P)) - X-o(P) - 0, CiXi(0"(P)) - Xi{P) = 0, CiX2(0"(P)) - X2{P) - 0, 

C2j/o(</'"(/'))-yo(/') = 0, C2yi(0"(P))-2/i(P)=O, and C2y2(0"(P)) - y2(P) = 0, 

for some nonzero constants ci,C2 G K. We know that there are only finitely many integers 1 < d < n 
and each value ciXi{4)'^{P)) — Xi{P) and C2j/i (P)) — yi{P) for < i < 2 can be factored into finitely 
many primes in K, so there can only be finitely many primes such that (j>'^{P) = P for d < n. □ 

The algorithm proceeds as follows to find a Q-rational periodic point with specified primitive period 
n on a given surface 5*: 

(1) Chose a prime p G Q. 

(2) Find all of the points in S{¥p). (Discussed in Section [3?T|) 

(3) Find all of the cycles in S{¥p). 

Exit if S is degenerate. 

(4) Save all points with m \ n for p. 

Exit if there are none. 

(5) Repeat for the next prime until the specified number of primes is completed. 

(6) Use the Chinese Remainder Theorem on the coordinates of the saved points to find P G 5'(Q). 
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There are two points worth further discussion: What ii S — V{L, Q) is singular and thus not a K3 
surface? Which points do you use for the CRT? 

First we address singularity. Since we are computing the forward image of every point in >5'(Fp) we 
can easily check for degenerate fibers since by [5] a fiber S* = p^^{a.) is degenerate if and only if 

G*{a) = G^(a) = G;{a) - H*,{a) = H*,{a) = H*,,_{ai) = 0. 

From [51 Proposition 2.5] this leaves the only possible singular points as fixed points of ctx and <7y. 
Therefore, if we are searching for a periodic point with primitive period n > 1 the singular fixed points 
have minimal effect on the algorithm. So we content ourselves to only check for degenerate fibers. When 
a surface 5* = V{L, Q) is found to have a periodic point we can then check to see if it is nonsingular. 

Determining which of the saved points to use in the Chinese Remainder Theorem is a more practical 
issue. We could try every possible combination of points with m | n modulo each prime p, however 
this may not be the most efficient strategy. Since we are only interested in finding a periodic point and 
not necessarily all periodic points and are willing to accept missing some periodic points for a faster 
algorithm, we choose only points from primes where there is a single point with m — n or primes where 
there is only a single point with m \ n. 

The above algorithm was used to find Wehler K3 surfaces with Q-rational periodic points of primitive 
period 1, . . . , 16. Table [T] and Table [3] give examples of surfaces with a point of each primitive period. 
Note that the surface with a 6-periodic point also has a 2-periodic point. 

Upon examining these tables, the reader will quickly notice that nearly all of the surfaces have coef- 
ficients contained completely in { — 1,0, 1}. Upon examining data from the first 30 primes for randomly 
generated 5* = V{L, Q), surfaces with larger coefficients had comparatively few periodic points: 48, 871 
out of 49, 200 surfaces (99%) with coefficients in [—128, 128] did not have a periodic point of any order. 
The remaining 329 surfaces may or may not have a periodic point but would require data from more 
primes to determine the existence of a periodic point. Whereas for surfaces with coefficients in { — 1, 0, 1}, 
7425 out of 19424 surfaces (37%) did not have periodic points of any order, with the remaining 11999 
surfaces requiring more data to determine the existence of periodic points. 

3.1. Finding p^^{P). The algorithm requires finding S'(Fp) efficiently. We do so by finding the 0, 1, 
or 2 points in p^^(x) for x G P^(Fp). We know from Call and Silverman [5l Corollary 1.5] that the two 
solutions {y,y'} satisfy the following equations: 





yoy'o 


= G§(x) 








yiy'i 


= GT(x) 








y2j/2 


= G^(x) 






yoy'i - 




= -ff?i(x) 






yoy'2 - 


^y'ov^ 








yiy'2 - 




= -^f2(x) 











= Gf(x)y^ 


-Hf2{y^)yiy2 - 


^GU^)yl 







= Gg(x)y2^ 




^G?(x)2/2 







= Gg(x)j/2N 


-^02(x)2/o2/2 - 





Solving this system is merely a matter of case by case analysis, so we leave the interested reader to 
examine the PARI/gp script ordersearch.gp or see [3 Section 6.2]. Note that the method is valid only 
for char(i^) = p ^ 2 due to using the quadratic formula for degree 2 polynomials, however, for F2 it is 
a simple matter to examine all possible points in P^(F2) x P^(F2). 
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n 


Point P and surface 5* = V{L,Q) 


1 


P = WO 11 \1 111 

L : [-1,1,1,0,-1,-2,0,1,0] 

Q ■■ [-2, -1, -1, 1, -2, 0, 1, -2, -1, 0, 1, -1, 1, 1, 1, 0, -2, -2, -1, -2, -1, -2, -2, 1, 
0,-2,-1,1,-1,-1,1,1,-2,1,1,1] 


2 


P — ff3 1 31 fl Oil 

L : [1,-3,2,0,1,-3,-1,-3,-1] 

Q : [-2, -4, -2, -2, -3, 2, 0, -2, 0, 1, -2, 0, 1, 2, 1, -4, 0, 0, 3, 2, -3, 1, -2, 1, -4, -2, 
1,1,2,2,2,-1,1,2,-4,0] 


3 


p — \\^ —1 01 \^ oil 

^ — [[i, i, UJ, [1, U, UJJ 

L : [0,-2,0,0,-2,1,1,-1,-2] 

Q : [0, 1, 0, 1, -1, -2, -1, 0, 0, 0, -1, -2, -2, 1, -1, -1, -1, 1, -1, -2, -1, -1, 0, 1, 
1, 0, -2, -2, 1, 1, -2, -2, -2, -2, 1, -2] 


4 


p — [[1 -il [1 _i nil 

L : [-1,-1,0,-1,0,-1,0,0,-1] 

Q : [-1,0,-1,0,0,-1,-1,0,-1,0,0,-1,-1,0,0,0,0,0,0,0,0,-1,0,-1,0,-1, 
0,0,-1,0,0,-1,0,-1,0,0] 


5 


P — ffl -1 01 fl 111 

L : [0,0,-1,-1,0,0,0,-1,-1] 

Q : [-1,-1,0,0,0,0,-1,0,0,0,-1,-1,0,0,-1,0,0,0,0,-1,0,-1,0,-1,-1,0, 
-1,0,-1,-1,-1,-1,0,0,0,-1] 


6 


Pfi = 111 OHl Oil A. = WO 1 OHl -111 

^ D LL 1 ^7 ^1 7 L 7 ^' "11 ^ LL"l 7 "1 7 L 7 "5 11 

L : [0,0,-1,-1,0,-1,0,-1,0] 

Q : [0,0,0,-1,0,-1,-1,0,0,0,-1,-1,-1,-1,0,0,-1,-1,0,-1,0,0,0,0,0,-1, 
0,0,-1,0,0,0,-1,0,0,0] 


7 


P= [[1,-2,1], [1,-1,-1]] 

L : [0,-1,0,-1,0,-1,-1,0,0] 

Q : [-1,-1,-1,-1,0,0,0,0,-1,-1,-1,-1,-1,0,0,0,0,0,0,0,-1,-1,0,0,0, 
-1,-1,-1,0,-1,0,0,0,-1,-1,-1] 


8 


P= [[1,0,0], [1,0,0]] 

L: [0,0,-1,-1,-1,0,0,-1,-1] 

Q : [0,-1,0,-1,-1,0,0,0,0,-1,-1,-1,-1,-1,0,0,-1,-1,0,0,-1,-1,-1,0, 
0,0, -1, -1, -1,0,0, -1, -1,0,0, -1] 



4. RiEMANN Zeta Function of a Wehler K3 surface 

Let p G Z be a prime and let = ^S{¥pm) be the number of Fp.Ti -rational points on S. Using 
Section [3.11 we can quickly determine the number of points in a fiber, and hence can efficiently determine 
Nm ■ We can then determine the Riemann zeta function of a surface S and give an example of a family 
of Wehler K3 surfaces that have Picard number 2. 

Definition 3. The Riemann zeta function of S is 

m=l / 



where exp denotes exponentiation. 
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Table 2. Primitive Periods 9-16 



n 


Point P and surface S ~ V{L, Q) 


9 


P ^ WO 1 01 \l Oil 

-t [^U, -L, VJJ, [-L, VJ, UJJ 

L : [-1,0,0,0,0,-1,0,-1,0] 

Q : [-1,0,0,-1,-1,-1,-1,0,-1,-1,-1,-1,0,0,-1,-1,-1,-1,0,-1,-1,0, 
0,0,-1,-1,0,-1,-1,0,0,-1,-1,-1,0,0] 


10 


P — ffl —1 01 fO 1 -111 

L : [0,-1,0,-1,-1,0,0,-1,-1] 

Q : [-1,-1,0,0,0,0,0,0,0,0,-1,-1,-1,0,0,-1,-1,-1,0,0,0,-1,-1,0,-1, 
0,-1,0,-1,0,-1,-1,-1,0,0,0] 


11 


p — \\^ n —11 f— 2 1 oil 

L : [-1,0,0,-1,-1,-1,-1,0,-1] 

Q : [0,-1,0,-1,-1,0,0,0,-1,-1,0,-1,0,0,0,0,-1,0,0,0,0,-1,-1,0,-1,-1, 
-1,0,0,-1,0,0,-1,-1,-1,-1] 


12 


P — \\2 1 -21 fO —2 111 

L : [-1,-1,-1,-1,-1,0,0,-1,0] 

Q : [0,-1,0,0,0,0,-1,0,0,-1,0,0,0,0,-1,-1,-1,-1,-1,-1,0,-1,-1,0, 
-1,-1,0,-1,0,-1,-1,0,-1,-1,-1,-1] 


13 


p — ffo 1 01 \^ -1 -111 

L : [0,-1,0,-1,-1,0,-1,0,-1] 

Q : [-1,0,-1,-1,0,0,-1,0,-1,0,0,-1,0,-1,0,0,-1,-1,-1,0,-1,0,0,0,0, 
-1,0,-1,-1,-1,0,-1,0,0,0,-1] 


14 


P = \\Q -2 -11 f-1 -3 211 
L : [0,-1,-1,-1,0,0,0,0,-1] 

Q : [0,0,-1,-1,0,0,0,-1,0,0,0,0,-1,-1,-1,-1,-1,0,-1,-1,-1,0,-1, 
0,0,-1,0,0,0,-1,-1,0,-1,-1,-1,0] 


15 


P= [[1,1,0], [1,-1,-1]] 

L : [-1,-1,0,-1,0,-1,0,-1,-1] 

Q : [0,-1,-1,0,0,0,-1,0,0,0,-1,0,0,-1,0,0,-1,-1,-1,-1,-1,0,0,-1,0, 
0,-1,-1,0,0,-1,-1,0,0,-1,0] 


16 


P= [[1,0,-1], [1,-1,0]] 

L : [-1,-1,0,-1,0,0,0,0,-1] 

Q : [0,0,-1,-1,0,0,-1,0,-1,0,-1,-1,-1,0,0,0,-1,0,-1,0,-1,0,0,-1,0, 
0,-1,-1,0,0,-1,-1,-1,0,0,0] 



The Riemann zeta function encodes much information about the set of points of S. Knowing the 
Ricmann zeta function therefore leads to knowing information about S. In particular, we produce a K3 
surface with Picard number two. One way to compute tlic Ricmann zeta function is to determine Nm for 
enough m. How many is enough depends on the variety. Since 5 is a K3 surface, it has Betti numbers 
&o = &4 = 1, ^1 = &3 = 0, and &2 = 22, consequently, we will sec that we need only the first 11 to 
determine the Riemann zeta function. 



4.1. The Riemann Hypothesis. The Riemann zeta function satisfies the Riemann Hypothesis, as 
shown by Dclignc 0. We have that dim(S') — 2 and that K ~ We are using p — i since the 
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Z{S, T) = exp 




Po{T)P2{T)Pi{T) 



algorithm from Section [Ql is not applicable in the case p = 2. Since 5 is a K3 surface, we have that 

m 

^ (i-r)(n-'i(i-«.T^))(i-9T)' 

Our goal is to calculate the 22 a^s. The method we use is to 

(1) Compute the NmS programmatically. 

(2) Write equation ([T]) as a power series in T. 

(3) Equate coefficients of the powers of T in the expansion of equation (m with the iV^s. 

From the above information, we arc able to calculate the power sums of the a^, from which we can 
compute the symmetric polynomials in the a^. Consequently, we can determine the polynomial with the 
ai as its roots. 

Fortunately, the 22 a^s are not all independent, and we show next that it is sufficient to determine 
only the first 11 NmS. The functional equation of the Riemann zeta function says that we can group the 
aiS into pairs which satisfy ata'^ = 9. Hence, we get 

22 11 11 

P2{T) = - a,T) = J|(l - a,r)(l - a[T) = J|(l - a,T + QT^), 

i—l i—1 i—1 

where 

9 , 

tti = a, H and aiOj = 9. 



Taking the natural log of both sides and expanding we get the following equations, where Pj = ^ 



11 3 

.-1 ai 



Pi 


= Ni- 


10 










P2 


= ^2 + 


- 116 










P3 


= ^3 + 


- 27Pi - 


730 








Pi 


= A^4 4 


36P2 - 


8344 








P5 


= A^5 4 


- 45P3 - 


405Pi - 


59050 






Pe 


= Ne^ 


54P4 - 


729P2 - 


515404 






Pi 


= Nr^ 


63P5 - 


II34P3 - 


F 5103Pi - 


4782970 




Ps 


= Ns^ 


72P6 - 


I62OP4 


f 11664P2 - 


- 43191064 




P9 




- 8IP7 - 


2187P5 - 


F 2I87OP3 - 


- 59049Pi - 


387420490 


PlQ 




f 90P8 - 


- 2835P6 


+ 36450P4 


- 164025P2 


- 3485485324 


Pll 




f 99P9 - 


- 3564P7 


+ 56133P5 


- 360855P3 


+ 649539Pi - 31381059610 



From these 11 power sums, it is possible to recover the first 11 symmetric polynomials Ci, . . . , Cn in the 
QiS using the Newton-Girard formulas [12], which relates the Pk to the Ci. We now have the polynomial 
whose roots are the a^, and we know that 



a^ — aiai + 9 = 0. 



Therefore, we can solve for 

± ^af - 36 



a. 



2 

So we determine the ai and can then determine the a,;. Therefore, we only need to determine the first 
11 N,nS to determine Z{S,T). 
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Table 3. Point Counting Data 



m 




1 


13 


2 


97 


3 


784 


4 


6877 


5 


60238 


6 


533440 


7 


4782322 


8 


43047613 


9 


387464230 


10 


3486563272 


11 


31382110828 



Remark. We could compute the Ci for p ^ i \i wc could feasibly determine the N,nS in this case. 

4.2. An Example Computation of the NmS. Thanks to the computation cluster at Amherst College 
for providing CPU time, we determine the first 11 for the Wehler K3 surface 

L ^ [1,0,0,0,1,0,0,0,1] 

Q = [2,2,2,0,1,1,1,2,1,1,2,1,0,0,2,2,2,2,0,1,0,2,0,2,2,0,2,1,1,1,0,1,0,0,0,0] 

The N,nS are shown in Table |3l From the N,nS we compute P2{T) to be 

PiiT) = (9T^ -6T+ 1) (3486784401^2° + 1162261467ri^ + 258280326T1® + 4304672iri^ 

- 143489077^^ - 143489077^^ - 4782969ri'' - bSlUlT^^ + 177U7T^^ + 118098r" 

+ 78732ri" + 13122T^ + 2187r® - 729T^ - 7297^ - - 27T^ + 9T^ + 6T^ + 3T + 1). 

Examining the we know from (TH Corollary 2.3] that the number of ^ which are roots of unity is 
an upper bound on the Picard number. In this example, the upper bound is then 2. However, we also 
know that the Picard number is at least 2, see [15]. Hence, the surface has Picard number exactly 2. 
Moreover, we have an infinite family of Wehler K3 surfaces with Picard number 2, since the Riemann 
zeta function over F3 is unchanged by adding multiples of 3 to the coefficients. 
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